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Impulse and sampled-data optimal control of heat equations, and 

error estimates 


Emmanuel Trelat* Lijuan Wangl Yubiao Zhangl 


Abstract 

We consider the optimal control problem of minimizing some quadratic functional over all possible 
solutions of an internally controlled multi-dimensional heat equation with a periodic terminal state 
constraint. This problem has a unique optimal solution, which can be characterized by an optimality 
system derived from the Pontryagin maximum principle. We define two approximations of this 
optimal control problem. The first one is an impulse approximation, and consists of considering a 
system of linear heat equations with impulse control. The second one is obtained by the sample-and- 
hold procedure applied to the control, resulting into a sampled-data approximation of the controlled 
heat equation. We prove that both problems have a unique optimal solution, and we establish precise 
error estimates for the optimal controls and optimal states of the initial problem with respect to its 
impulse and sampled-data approximations. 
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1 Introduction and main results 

1.1 The context 

There is a vast literature on numerical approximations of optimal control problems settled for parabolic 
differential equations. In the linear quadratic regulator (LQR) problem, many results do exist concerning 
space semi-discretizations of the Riccati procedure. We refer to [1, 2, 14, 16, 27, 31] for general results 
showing convergence of the approximations of the Riccati operator, under assumptions of uniform expo¬ 
nential stabilizability, and of uniform boundedness of the sequence of approximate Riccati solutions. In 
[1, 22, 27], these sufficient conditions (and thus, the convergence result) are proved to hold true in the 
general parabolic case and for unbounded control operators. Note that, in such LQR problems, the final 
point is not fixed. When there is a terminal constraint the situation is more intricate, because things may 
go badly when discretizing optimal control problems in infinite dimension, due to interferences with the 
mesh that may cause the divergence of the optimization procedure when the mesh size is going to zero. 
These interferences are stronger when the terminal constraint has infinite codimension, in spite of strong 
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dissipativity properties of parabolic equations. For the optimal control problem of minimizing the 
norm of the control (corresponding to the celebrated “Hilbert Uniqueness Method”), one can find results 
on uniform exact controllability and/or observability of discretized control systems in [4, 5, 7, 8, 18, 38] 
(see also references therein), for different discretization processes on different parabolic models. It can 
be noted that uniformity requires in general to add some appropriate viscosity terms in the numerical 
scheme. Besides, when the convergence is ensured, it is important to be able to derive error estimates 
which are as sharp as possible, and we refer the reader to [10, 12, 13, 15, 20, 21, 26, 32, 33, 34] for 
situations where Galerkin finite element approximations are used. 

In many cases impulse control is an interesting alternative, not only to usual discretization schemes, 
but also in order to deal with systems that cannot be acted on by means of continuous control inputs, as it 
often occurs in applications. For example, relevant controls for acting on a population of bacteria should 
be impulsive, so that the density of the bactericide may change instantaneously; indeed continuous control 
would enhance drug resistance of bacteria. For more discussions and examples about impulse control or 
impulse control problems in infinite dimension, we refer the readers to [3, 36, 35] and references therein. 
It is also interesting to note that impulse control is as well an alternative to the well known concept of 
digital control, or sampled-data control, which is much used in the engineering community. 

To the best of our knowledge, error estimates for impulse approximations or for sampled-data approx¬ 
imations of an optimal control problem settled with partial differential equations and with continuous 
control inputs have not been investigated. 

In this paper, we consider the problem of deriving precise error estimates for impulse approxima¬ 
tions and for sampled-data approximations of a linear quadratic optimal control problem settled for an 
internally controlled linear homogeneous heat equation with periodic terminal state constraint. The lat¬ 
ter periodicity requirement is motivated by the fact that steady solutions and periodic solutions are of 
particular interest when considering parabolic differential equations. 

1.2 Definition of the optimal control problems 

Let ^ 1 be an integer, let H C be a bounded open set having a boundary dfl, let w C fl be an 
open non-empty subset, and let T > 0 and yd S L^(0, T; L^(H)) be arbitrary. Throughout the paper, the 
norm in L^(H) is denoted by jj jj. 

The optimal control problem (OCP). We consider the optimal control problem (OCP) of mini¬ 
mizing the functional 

= \\y - VdW^ At + ^ IlMpdt (1.1) 

over all (y, u) S L^(0, T; il^(H) fl Hq{Q)) fl il^(0,T; L^(H)) x L^(0, T; L^(H)) such that 

( dty-Ay = XujU in Hx(0,T), 

y = 0 on dnx{0,T), (1.2) 

i 2/(0) = y(T) in n. 

Here, Xoj designates the characteristic function of w, and y and u are functions of (t, a;). 

We have the following facts: 

• Given any u £ L^(0,T;L^(r2)), there exists a unique solution y £ L^(0,T; 71^(H) fl HQ{n)) fl 
H^{0,T-L^{n)) of (1.2). 

• The problem (OCP) has a unique optimal solution {y*,u*). 
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• By the Pontryagin maximum principle (in short, PMP; see [24]), which is here a necessary and 
sufhcient condition for optimality because the problem is linear quadratic, this minimizer is char¬ 
acterized by the existence of p* G i/^(0,T;fl L^(0,T;niJQ(r2)) such that 

dty* - Ay* = XujU* in r2x(0,r), 

y*=0 on dnx{0,T), (1.3) 

y*{0)=y*{T) in fl, 

dtp* + Ap* = y* - yd in flx(0,r), 
p*=0 on anx(0,T), 

p*{Q)=p*{T) in fl, 

and 

u* = Xu:P* in n X (0,T). 

These three claims are easy to establish, but for completeness they are proved in Section 2.1. 

We are next going to design an approximating impulse optimal control problem (lOCP)n, and an 
approximating sampled-data optimal control problem (SOCP)n, for a linear heat equation with periodic 
terminal state constraint. Both problems have as well a unique solution, to which we will apply the 
PMP. We will then establish error estimates between the optimal solutions of (OCP) and, respectively, 
(IOCP)„ and (SOCP)„. 


(1.4) 

(1.5) 


The approximating impulse optimal control problem (IOCP)„. Let us define the approximating 
impulse optimal control problem (lOCP)n, for n ^ 2. We set 

hn = T/n, Ti=ihn, z = 0, l,...,n, 

n 

X = = z = l,2,...,n, 

2=1 


and we define the functional Jn '■ X x (L^(ri))” ^ —>■ [0, +oo) by 

l|Wi-l,nf ^ , (1-6) 

for w = (2/1,n, 2/2,n, • ■ • , 2 /n.n) G X and Un = (Mi,„,it2,n, ■ • .,Un-i,n) G (^^(fl))"”^ Here and throughout, 
II • II designates the norm in L^(H). Accordingly, the inner product is denoted by (•, •). 

We consider the impulse optimal control problem (lOCP)n, consisting of minimizing the functional 
Jn over all possible {Yn,Un} G AT x (L^(H))”“^ such that 

{ ^ty^,n - Ayi^n = 0 in 

yi,n = 0 on 

2/i,n(ri_i) = 2/*-i,n(ri_i) -I- Xi^u^-i,n in 
2/1, n(0) = yn,n{T) in 

Proposition 1. For every n ^ 2, the optimal control problem (lOCP)n has a unique solution (Y*,U*), 
with Y* = {yln,yln^---,yn,n) U* = uWi„). The Optimal solution (Y*,U*) of 

(IOCP)„ is characterized by the existence of p^ G L^(0, T; iL^(fl) niLo(f2)) fl iL^(0, T; L^(H)) such that 


H X 1 ^ ^ n, 

an X (ri_i,Ti), 1 ^i^n, 
H, 2 ^ ^ n, 

H. 


Jn{Yn,Un) = 


= i w 


•'n-l 


|2/*,n - l/dll^dA 


■E 

2=2 


dty*,n - ^Vln = 0 

in 

H X (ri_i,ri), 

1 ^ ^ n, 

Vlu = 0 

on 

an X {Ti-i,Ti), 

1 ^ ^ n 

Vlniji-l) = y*_i^n{n-l) + 

in 

FL, 2 ^ i ^ n, 


2 /l*,n( 0 ) = y*njT) 

in 

H, 



( 1 . 8 ) 
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( dtPn + Ap* =yl-yd in Fix (0,T), 
p; = 0 on do. X (0,T), 

Ik(0)=p);(A) in fl, 

(1.9) 

and 

Wi-l.n = KXi^Pniji-l), 2 4 7 4 77 , 

(1.10) 

with 

p:(0) = p:(T), 

(1.11) 

where y* G L°^{0,T; L^{n)) 

is defined by 



yUt) = ylnit), te{n-i,n], 1 4 f 4 n. 

(1.12) 


and G C{[n_i,n]; L'^{n)), 1 ^ i < n. 

Proposition 1 is proved in Section 2.2. 

Remark 1. We could also consider the corresponding impulse version of the optimality system (1.3)- 
(1.4)-(1.5), but then its well-posedness would be hard to prove, and therefore, obtaining error estimates 
in such a way seems difficult. 


The approximating sampled-data optimal control problem (SOCP)^. Let us now define the 
approximating sampled-data optimal control problem (SOCP)n, for n ^ 2, by performing the usual 
sample-and-hold procedure on the control function. This consists of freezing the value of the control over 
a certain horizon of time, usually called sampling time. In other words, we replace the control function 
u G L^(0, T; with a control that is piecewise constant in time, with values in We keep 

the same notations as in the definition of (IOCP)„, and we assume that the sampling time is equal to 
hn = T/n. Recall that we have set Ti = i hn, for f = 0,..., n. We consider the class of sampled-data 
controls /„ G L^(0, T; L^(n)) defined by 

fn(t)=Vi^n, Vt G (Ti_i,Ti], 1 < I ^ 71, (1.13) 

where Vi^n G L‘^{Vl) for every i G {1,... ,n}. This class of controls is therefore identified with (L^(r2))". 

Recall that the functional J is defined by (1.1). We consider the sampled-data optimal control problem 
(SOCP)n, consisting of minimizing the functional 


over all (2/„,14) G L^(0, T; iJ^(fl) niL,)(n)) niL^(0, T; L^(n)) x with 14 = (ui,„,... ,u„,„), such 

that 

dtVn-^yn = Xu>fn in f2x(0,T), 

Un = 0 on dnx(0,T), 

2/«(0) = yn(T) in n, 

where /„ G L^(0, T; L^(r2)) is the sampled control defined by (1.13). 

Proposition 2. For every n ^ 2, the optimal control problem (SOCP)^ has a unique solution {y*, V*), 
with V* = (u* ..., u* „). The optimal solution o/(SOCP)„ is characterized by the existence 

of PI G niJi(0,T;L2(f7)) such that 

dtvl - ^y*u = Xujfn n X (0, T), 

y: = 0 077 dflx(0,T), (1.14) 

yn(.o)=y*{T) in fl, 


J{yn,fn) = ^(^j^ ||7/„-7/df dt- 


^Ei 
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and 


f dtpl 

+ Ap; = yl-yri 

in 

Dx (0,r), 


\ P*n = 

0 

on 

dXl X 

(0,T), 

(1.15) 

1 p);(o) 

= Pn{T) 

in 

D, 



* 

= -^Xu^ I Pnit)' 

J Ti-\ 

dt, 

1 ^ i 

< n. 

(1.16) 

the (optimal) sampled-data control given by 


mt) 

= Kn, Vt G (Ti-i 

,Ti], 

1 ^ i 

^ n. 

(1.17) 


Since the proof of Proposition 2 is similar to the one of Proposition 1, we do not provide any proof in 
the present paper. It is interesting to note that the optimal sampled-data control /*, defined by (1.17), 
is given by time-averages of the adjoint state p* over the time-subdivision defined by the sampling time 
hn (see (1.16)). This fact has been proved as well in a more general context in [6]. 

Remark 2. It is clear that the sampled-data optimal control problem (SOCP)„ may be considered as 
an approximate version of (OCP), but it is less clear, at least intuitively, for the impulse optimal control 
problem (IOCP)„. Before establishing precise error estimates in the next section, let us provide a first 
intuitive explanation. Firstly, any continuously distributed control may be discretized by the sample-and- 
hold procedure, leading to the sampled-data control with vt € L‘^{Vl). Secondly, 

this sampled-data control can be seen as an approximation of an impulsive control, in the sense that 

1 . 

T-X(0,hn) -t <5{t=0} 

'^n 

in the distributional sense, where is the Dirac mass at t = 0 (note that we have as well the 

convergence of the corresponding solutions, see Lemma 2 further). Denoting by y{u) the solution of 
(1.2), we have, noting that Ui-i^n — hnu{ihn), for f = 2,..., n. 


y{u) ~ y 



)] - y 


\i^2 


h„) ~ y 


/ . '^{‘= 6 - 1 ) hr 
^i=2 


}^2—l,n 


and 

n 1 ^ 

ll^lli2(0,T;L2(O)) — y]] hn\\u{i h„)||^ ~ ||ui_i,„||^. 

i=l " i=2 

Here, y (X]r =2 '^{*=(*-1) i® corresponding impulsive solution, and we have moreover J{y, u) ~ 

JniYn,Un). 


1.3 Error estimates 

We keep all notations introduced in Section 1.2. The main results of the paper are the following. 

Error estimates for the impulse approximation. 

Theorem 1. We set 


Kit) = 


t G (ri_i,ri], 1 ^ f ^ n. 




= 0 . 


(1.18) 
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Then there exists G{T) > 0 such that 


11^* “ '*^nllL2(0,T;L2(n)) ^ G{T)hl/^\\yd\\ L^(o,T-,L^(Q)), (1-19) 

and 

|J„(y;,C7:)-J(j/*,w*)KC(T)/iy2||yd|li.(o.T;L^(a))- (1-20) 

For every p G [2, +oo), there exists G(T,p) such that 

\\y* - yn\\LP{o,T-L^(n)) ^ G{T,p)hl/^\\yd\\L^(o,T-,L^(n)), pS [2,+oo). (1-21) 


Moreover, the constant G{T) and G(T,p) are independent of n and of yd. 


Theorem 1 is proved in Section 2.3. Note that we have assumed that p < +oo in the statement. If 
p = +00 then the situation is more complicated, and we have the following result. 


Theorem 2. ITe assume that the subset lo of has a G^ boundary. Let q G (l,+oo) be arbitrary. If 
Lo ^ Ft, then 


f G{T)hn'^^\\yd\\L'^(Q^T-L'^{Q.)) for N > 3, 

l|y* - 2 /JlllL°°( 0 ,T;L 2 (n)) >1 G{T,q)hl!‘^‘^\\yd\\L'^(o^T-L'^(n)) for N = 2, 

y C'(T)ft,y"^||?/d||i 2 (o,T;L 2 (o)) for N = l. 


( 1 . 22 ) 


If uj = Ft, then 

\\y* - 2/))IU“(0,T;L2(n)) < G{T)hl/‘^\\yd\\L^{o,T-L^{Q))- 
The constant G(T) and G(T,q) are independent of n and of yd. 

Theorem 2 is proved in Section 2.4. 

Remark 3. The above error estimates are much easier to obtain when the control domain oj is equal to 
the whole domain Ft, that is, when uj = Ft. But in this case the optimal control problems (OCP) and 
(lOCP)n have little interest. Actually, the main difficulty in obtaining our results is due to the fact that, 
if w C 0, then the function Xu: is not smooth and the function Xu>Pn{Ti-i) iii (1-10) is not in Rq)!!). In 
the proofs of Theorems 1 and 2 (see Section 2), in addition to more or less standard functional analysis 
arguments, to overcome the abovementioned difficulty, we use smooth regularizations of the characteristic 
function x^ji gradient of which we have to estimate in a refined way in some appropriate norm. 
Of course, this gradient blows up as the regularization parameter tends to zero, but fortunately there is 
some room to design appropriate regularizations, with adequate blow-up exponents (which we compute 
in a sharp way) that can be compensated elsewhere in the estimates, using Sobolev embeddings and 
usual functional inequalities. Using this approach, and deriving nonstandard estimates for the linear 
heat equation, we ultimately establish the desired error estimates. 

Remark 4. In Theorem 2, if w = fl (trivial case, according to Remark 3) then the order of convergence 
of the state is 1/2, and we conjecture that it is sharp. ^ If w C O, then we have obtained the error estimate 
(1.22) but we conjecture that it is not sharp, and that the order of convergence 1/2 should hold true as 
well. 

^Actually, we are able to prove that the exponent 1/2 is sharp in the estimate given in Lemma 2 (in Section 2.3) and in 
Lemma 5 (in Section 2.4.3), in the case where uj = Q and p = 2. We do not provide the proofs of these facts here. 
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Error estimates for the sampled-data approximation. 
Theorem 3. There exists C{T) > 0 such that 


Ik* - /nllL=(0.T;L2(n)) 



1/2 


U - Vi 


'dt ^ C'(T)/i„||yd||i2(o_'r;L2(n)), 


(1.23) 


Ik* “ J/nllc([0,T];ffi(f2)) + Ik* “ VnW L^{ 0 ,T-m (Q,T-L2{Q)) ^ C!{T)hn\\yd\\L^{ 0 ,T-,L^(n)), (1-24) 

and 

\Jiy*, u*) - Jifn. f:)\ ^ CiT)hJyd\\h^o.T-,mn)y ( 1 - 25 ) 

Moreover, the constant C{T) is independent of n and of yd- 

Theorem 3 is proved in Section 2.5. 

Remark 5. It is interesting to note that the error estimates are better for the sampled-data approxi¬ 
mation than for the impulse approximation. For instance, the control error estimate is of order 1 for the 
sampled-data approximation, but is of order 1/2 for the impulse approximation (in norm). This is 
not surprising because, as explained in Remark 2, the sampled-data optimal control problem (SOCP)n 
can easily be recast as a classical approximation of (OCP), in the sense that the class of admissible 
sampled controls is a subset of the class of admissible controls of (OCP). In this sense, obtaining the 
error estimates of Theorem 3 could be expected. In contrast, the set of unknowns (T„, [/„) for the impulse 
optimal control problem (lOCP)n is not a subset of the set of unknowns {y, u) for (OCP). This explains 
why the derivation of error estimates for (lOCP)n is much more difficult. 


1.4 Further comments 

We have established error estimates for the optimal controls and states of impulse approximations and 
of sampled-data approximations of a linear quadratic optimal control problem for a linear heat equation, 
with internal control, and with periodic terminal constraints. To our knowledge, this is the first result 
providing such convergence results and estimates, in an infinite-dimensional context. Many questions are 
open, that are in order. 


Terminal constraints. Here, we have considered periodic terminal constraints. This condition is 
instrumental in order to obtain existence and uniqueness results and to be able to derive a PMP (see, in 
particular. Lemma 1 in Section 2.1). But it is of course of interest to consider other terminal conditions. 
For instance, one may want to consider the problem (OCP) with the fixed terminal conditions y{0) = y^ € 
Lf{Tl) and y{T) = 0. It is well known that this exact null controllability problem admits some solutions, 
without any specific requirement on the (open) domain of control w. But it is well known too that the 
final adjoint state coming from the PMP lives in a very big space that is larger than any distribution 
space. This raises an important difficulty from the functional analysis point of view, preventing us from 
extending our analysis to this setting. 

Moreover, when considering more general equations (see the next item), if one considers an infinite- 
codimensional state constraint then it is well known that the PMP may fail (see [24]), and then in this 
case even the basic fact of establishing an optimality system may raise some impassable obstacles. 
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More general evolution equations. We have considered the linear homogeneous heat equations. The 
question is open to extend our analysis to more general parabolic equations, of the kind dty = Ay + XujU, 
with A : D(A) —^ L^(r2) generating an analytic semi-group. For instance, one may want to replace 
the Dirichlet Laplacian with a general elliptic second-order differential operator, with various possible 
boundary conditions (Dirichlet, Neumann, Robin, mixed), or with the Stokes operator. It is likely that 
our results may be extended to this situation, but note anyway that, in our proofs, we repeatedly use the 
fact that we have Dirichlet conditions. 

More generally, the question is open to investigate semilinear parabolic equations, of the kind dty = 

+ f{y) + XtjM. Even when A is the Dirichlet Laplacian, this extension seems to be challenging. 

The case of hyperbolic equations is another completely open issue. Certainly, the first case to be 
investigated is the wave equation: in that case one replaces the first equation in (1.2) with the internally 
controlled linear homogeneous wave equation dtty = Ay -|- In this case, it is well known that exact 
controllability holds true under the so-called Geometric Control Condition on {u!,T). What happens for 
the corresponding impulse model is open, and is far from being clear (see [17] for results in that direction). 
Also, many estimates, that are quite standard for heat-like equations and that we use in this paper, are 
not valid anymore in the hyperbolic context. 

More general control operators. In this paper, we have consider an internal control. Writing the 
control system in the abstract form dty = Ay + Bu, this corresponds to considering the control operator 
B defined by Bu = XcjM. In this case, the control operator is bounded, and we implicitly use this fact 
in many places in our proofs. We expect that our results can be extended to more general classes of 
bounded control operators, but the case of unbounded control operators seems much more challenging. 
For instance, what happens when considering a Dirichlet boundary control is open. 

Time-varying control domains and optimal design. Another open question is to derive our error 
estimates for time-varying control domains. Note that control issues for wave equations with time-varying 
domains have been investigated in [23]. In our context, this means that we consider a control domain 
uj{t) depending on t in (1.2). In this case, the definition of the approximating impulse control system 
(1.7) must be adapted as well, by considering uj{Ti-i) at time Ti_i. It is likely that our main results may 
be at least extended to the case where w(t) depends continuously on t. The general case is open. 

Related to this issue, is the question of determining how to place and shape “optimally” the control 
domain. Of course, the optimization criterion has to be defined, and we refer to [28, 29, 30] where optimal 
design problems have been modeled and studied. In the context of the present paper, we could investigate 
the problem of designing the best possible control domain such that the constants appearing in our error 
estimates be minimal. 

Let us be more precise and let us define the open problem. Given any T > 0, and any open subset 
uj of D, Theorem 1 asserts that there exists a constant C'(T) > 0 such that the error estimates (1.18), 
(1.19) and (1.21) (with p = 2, for instance) are satisfied. Since this constant depends on w, we rather 
denote it by Ct{^)- Given a real number L G (0,1), we consider the optimal design problem 

inf C' 7 ’(a;), 

that is, the problem of finding, if it exists, the best possible control subset having a prescribed Lebesgue 
measure, such that the functional constant in the error estimates is as small as possible. This prescribed 
measure is L|fl|, that is, a fixed fraction of the total volume of the domain. We stress that the set of 
unknowns is the (very big) set of all possible measurable subsets of fl of measure L|D|. It does not share 
any good compactness properties that would be appropriate for deriving nice functional properties, and 


thus already the problem of the existence of an optimal set is far from obvious. However, following [30] 
where similar optimal design problems have been investigated in the parabolic setting, we conjecture that 
there exists a unique best control domain, in the sense given above. Proving this conjecture, and deriving 
characterizations of the optimal set, is an interesting open issue. 

Note that, less ambitiously than the problem above, one could already consider simpler optimal design 
problems, where the problem consists, for instance, of optimizing the placement of a control domain having 
a prescribed shape, such as a ball: in this case the set of unknowns is finite-dimensional (centers of the 
balls). 

Impulse Riccati theory. In the present paper, we have considered a problem within a finite horizon 
of time T. It would be interesting to consider the optimal control problems (OCP), (lOCP)n and 
(SOCP)n in infinite horizon, that is, when T = -|-oo. In this case, the optimal control solution of (OCP) 
is obtained by the well known Riccati theory (see, e.g., [37]), which gives, here, 

u = XLjE{y-yd), 

where E G L(L^(H)) is the unique negative definite solution of the algebraic Riccati equation 

AE + EA + ExujE = id. 

For the approximating impulse problem (lOCP)n, up to our knowledge, the Riccati procedure has 
not been investigated. In other words, up to now there does not seem to exist a Riccati theory for impulse 
linear quadratic optimal control problems in infinite dimension. Developing such a theory is already a 
challenge in itself. Assuming that such a theory has been established, the next challenge would be to 
establish as well the corresponding error estimates on the control and on the state, as done in our paper. 

For the approximating sampled-data problem (SOCP)„, few results exist in the literature. In [31] the 
authors have established a convergence result (which can certainly be improved, by combining it with the 
more recent results of [22, 27] for instance), but we are not aware of any result providing error estimates 
as in our paper. 

For impulse systems in particular, such a theory would certainly be very useful for many practical 
issues, because, as already mentioned, impulse control may be an interesting alternative to discretization 
approaches, or to sample-and-hold procedures, which is sometimes better suited to the context of the 
study. Notice that, although the theory of space semi-discretization of the Riccati procedure is complete 
in the parabolic case (but not in the hyperbolic case when the control operators are unbounded), to 
our knowledge the theory is far from complete for infinite-dimensional sampled-data control systems. 
Therefore, with respect to sample-and-hold procedures, this is one more motivation for developing an 
impulse Riccati theory and its approximations. 


2 Proofs 

2.1 Preliminaries 

Existence and nniqueness. We start with an easy existence and uniqueness result, together with 
a useful estimate. Throughout the paper, we denote by the semi-group generated by the 

Dirichlet-Laplacian on L^(H). 
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Lemma 1. Let T > 0 be arbitrary. Let f € L^(0,T; L^([})). Then the equation 

( dty — Ay = f in Ll x {0,T), 

J 2 / = 0 on dLtxiO,T), (2.1) 

i 2/(0) = y(T) in n 

has a unique solution y G H^{0,T; L'^{Ll)) D L^{0,T; H^{Ll) DHqITI)). Moreover, there exists C{T) > 0, 
not depending on f and on y, such that 

ll2/|lc([O.T];ffi(n)) + ll2/||ffi(0.T;L2(n))nL2(0.T;ff2(O)nffi(n)) ^ C(T)\\f\\i^ 2 (^Q^rp.j^ 2 (^Q^y (2.2) 

Proof. As a preliminary remark, we recall that, given j/q S L^{Ll) and / S Lf{Q,T-,Lf{Ll)) arbitrary, 
there exists a unique weak solution y{-;yo,f) S L^(0, T; iJQ(n)) D H^{0,T; H~^{n)) C C{[0,T]] L'^{Lt)) 
of dty — Ay = / in x (0, T), with y = 0 along dQ x (0, T), such that j/(0) = yo (see [25] for instance). 
Here, “weak” means that the differential equation is written in Moreover, if yo G Hq{LI), then 

2 /(-; 2 / 0 , /) e i"(0,T; H^in) n H^m n h\o,t- L\n)) c ^([o, T]-Hl{n)). 

Given any / S L^(0, T; L^(n)), let us prove the existence and uniqueness of a weak solution of (2.1). 
Since ||e^^||£(L 2 (Q)_i 2 (Q)) ^ < 1, it follows that (/-e^^)”^ exists and || (/-e^^)”^ ||£(L 2 (n)_i 2 (Q)) < 

(1 — where —Ai < 0 is the first eigenvalue of the Dirichlet Laplacian. Now we define 

yl = {I- e^^)-i r e(^-0^/(t) dt, (2.3) 

^0 

and 

2/(^;yo>/) = e*'^2/o + / e(‘"")^/(s)ds, te[0,r]. (2.4) 

Jo 

Then^Q G L^(f2) and ?/(•; ?/q ,/) G L^(0,T;i?g(fl))niJ^(0,T;i?“^(r2)) is the weak solution of c/jj/—Ay = / 
in 11 X (0,T), with y = 0 along dft x (0,T), such that y(0) = y^. Using (2.3) and (2.4), we have 

yiT;yLf) = e^^y^+ [ e^^-*^^f{t)dt = e^^{I-e^^)-^[ e(^-‘)^/(i)dt+ / e(^-‘)^/(i) dt 
^0 Jo Jo 

= {1- r e(^-‘)^/(t) dt = y(0; y^, /), 

^0 

which gives the periodicity requirement. Hence y(-;yo,/) is a weak solution of (2.1). 

Now, if yi and y 2 are two weak solutions of (2.1) associated with /, then 

^^ll 2 /i( 0 - 2 / 2 ( 0 lli 2 (n)+^|Vyi(t)-Vy 2 (t)pdx = 0, a.e. tG(0,T). 

Integrating the latter equality over (0, T), we deduce from the periodicity condition that yi = y 2 - There¬ 
fore the weak solution is unique. 

It remains to prove that the weak solution y of (2.1) actually belongs to L^{0,T; fl L[q{LI)) fl 

H^(0,T;L^(H)) and to prove the estimate (2.2). Using the preliminary remark, we have y(T) G Hq{LI), 
and since y(0) = y{T), it follows that y(0) G iJo(fl)- Therefore y G L^(0,T;i?^(H) fl iJg)!!)) fl 
H^{d,T\Lp'{Ll)). Now, multiplying the differential equation by 2y and integrating over 11, we get that 

^||yf-f2||Vyf = 2(/,y), a.e. lG(0,r). (2.5) 
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Using the Poincare inequality ||(/7|| ^ C||V(/3||, valid for any ip G combined with (2.5) and the 

Young inequality, we infer that 

^||yf+ 2||Vyf <||Vyf+ C||/f, a.e. t & {0,T). 

Here and throughout, C designates a generic positive constant only depending on H. Integrating over 
(0,T), we obtain that 


liVyf dt 


■dt. 


( 2 . 6 ) 


Besides, multiplying the first equation of (2.1) by —2tAy and integrating over H, we have 


t-||Vyf + 2t||Ayf = -2t{f,Ay) < t||Ayf + t||/f. 


Integrating again over (0,T), we obtain that 


T’l|Vy(T)f < TllVyrdt + T Tn/f dt, 
^0 ^0 


which, combined with (2.6) and the third equation of (2.1), gives ||Vy(0)|p ^ /o^ll/lPdt. This, 
together with the first and the second equations of (2.1), implies that 

ll2/llc([o,T];i/i(n)) + ll2/llffi(o.T;L2(n))nL2(o^T;ff2(n)nHi(n)) 

^ C'(||Vy(0)|| + ||/||l2(0,T;L2(O))) < ^ rp - ^ll/llL2(o.T;L2(n))- 


This completes the proof. 


□ 


Optimality system (PMP). The proof of existence and uniqueness of an optimal solution of (OCP) 
is easy. Since it is similar, but simpler, than the proof of Proposition 1 further, we skip it. 

Let {y*,u*) be the optimal solution of (OCP). For any v G L^(0, T; L^(f2)) and A G M \ {0}, let y\^v 
be the solution of 


Setting z 


{ dty\,v - Ayx,v = Xui{u* + Xv) 
y\,v = 0 

2/A,t,(0) = yx,v{T) 

, we have 


in H X (0,T), 
on dQ X (0, T), 
in n. 


dtz - Az = XujV 


z = 0 

z(0) = z(T) 


in H X (0,T), 
on dfi X (0, T), 
in fl. 


Moreover, by definition, J{y\,v,u* + Xv) > J{y*,u*), for every A ^ 0, and hence 



yd)zdx dt + 


*vdxdt = 0. 


(2.7) 
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Let p be the solution of 


dtp + Ap = y* - yd in 
p = 0 on 

p{0) = p(T) in 


This, together with (2.7), yields that 

0 = / [ {y* - yd)zdxdt 

Jo Jn 


u*vdxdt 


10 Jn 


n X (o,T), 
dn X (o,r), 
n. 


f I 


,p + M*)z;da;dt. 


Hence u* = Xi^P- This gives the PMP for the problem (OCP). 


2.2 Proof of Proposition 1 

Let us first prove the existence and uniqueness of a solution of (lOCP)n- According to the beginning 
of the proof of Lemma 1, we first recall that for Un = {ui^n, ■ ■ ■ ,Un-i,n) G we say that 

Yn = iyi,n, ■ ■ •, yn,n) G A is a weak solution of (1.7) if 


{dtyi,n{t),p)H~^{ 0 ),Hl(si)+ j^"^yi,n{t)-ypdx = d a.e. te(Tj_i,Ti), 1 < f ^ n, 

for each p G Ilglfl) (this means that the differential equation is written in iJ“^(H)) and yi,n(Ti-i) = 
yi-i,n(Ti-i) + XLjUi-i^n, for 2 ^ z ^ 71, and yi^nid) = yi^n{T). It is then a standard fact that (1.7) has a 
unique weak solution. 

Let d* = inf J„(y„, 17„) ^ 0, where the infimum is taken over all pairs {Yn,Un) G X x (L^(fl))”“^ 
satisfying (1.7). By definition, there exists a sequence (T„,m, 17„,m)m>i, with Yn^rn = (j/i,n,m)i<i^n and 
Un,m = (■Ui,n,m)i<i^n-1 Satisfying (1.7), such that 

d* ^ Ju{Yn,m, Un,m) ^ d* + -. ( 2 . 8 ) 

m 

Integrating the equations given by (1.7), we get 

yi,n,mit) = e*^yi^n,mi0), te[0,ri], 

i 

yi^n,m(d^ — ^ 2/l,n,m(0) T ^ ^ ^ ^ G [Tj_i,T 2 ], 2 ^ Z ^ 7Z, 

(2.9) 

n 

yi,n,m(0) — ^ yi,n,m(0) T ^ ^ —l,n,m. 

i=2 

Using (2.8) and the third equality of (2.9), we infer that 


||j/l.n.m(0)|KC, (2.10) 

for every m ^ 1.^ Here and throughout the proof, C designates a generic positive constant not depending 
of m. Multiplying the first equation of (1.7) (written for yi^n,m) by 2yi^n,m and integrating over H x 
(Ti_i,t), we obtain that 


|| 2 /i.n,m(i)r + 2 / ||Vyj,„,m(s)|pds < ||z/i,„,m(Ti_i)|p, Vt G [ri_i,ri], 1 ^ z < n, 

Jn-i 

^Here, the norm is used. For yi,n,m{0), we may wish to consider the Hq{Q) norm. But since yi,n,m{'^i—i) = 
yi-l,n, 7 n{n-l) ^Xuj'^i-I,n and (^) {2 ^ i ^ n), it follows that yi^n,min-l) € for 2 ^ ^ n. Hence the 

Hq{Q,) norm does not seem to be useful. 
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which implies that 

||yi,n,m llc([Ti_i,ri];L^(n)) + ||2/i,n,m||Xj ^ — l) II5 1 ^ ^ 7Z. 

This, together with (2.8), (2.10) and the third equation of (1.7), gives 

n n 

y^(||2/i,ra,m||c([Ti_i,Ti];L^(n)) + || 2/z,n,m || Xj ) + E ll'^z — l,n,7n|| ^ C. 

2—1 2—2 

Hence, up to some subsequence, we have 

yi,n.m^yln Weakly in Xi, strongly in r^; L^(ri)), 1 ^ i ^ n, 

and 

Wi-I,n.m Ui-i,n weakly in 2 ^ i n, 

for some Y* = (?/* „, y* , y* J £ X and U* = (u^ ..., G Passing to the 

limit in (2.8) and in (1.7), it is clear that {Y*,U*) is an optimal solution of (lOCP)n. 

The uniqueness follows from the strict convexity of the functional J„ : (L^(r2))"“^ —>■ [0, +oo) defined 
by JniUn) = JniYn,Un), where Yn is the unique solution of (1.7) corresponding to Un- 

Let us now prove the characterization of the optimal solution given in the proposition. 

We assume that {Y*,U*) is the optimal solution of (IOCP)„. Let us prove the existence of the adjoint 
state. The argument goes by perturbation of the optimal solution. Given any Un = (ui.n, U 2 ,n, ■ ■ ■, Un-i,n) S 
(L^(n))”“^ and any A G (0,1), we set 


Un.X = K+ XUn- 


( 2 . 11 ) 


Let Yn^x = (2/1, n, A, ?/2.n,A, • ■ •, J/n.n.v) be the solution of 


^yi,n^X — 0 

in 

D X (Tj_i,Ti), 

1 ^ Z ^ 7Z, 

yi,n,X — 0 

on 

5D X {Ti-i,Ti), 

1 ^ Z ^ 7Z, 

yi,n,x{'^i-l) = yz-l,n,A(Ti-l) + Xa;«-l,n + ><Ui-i^n) 

in 

D, 2 ^ z ^ 7Z, 


yi,n,A(^) ~ yn,n,A(T') 

in 

n. 



Setting Zi^n = 1 ^ i ^ n, we have 


in nx(Tj_i,rj), 
on dfl X {Ti-i,Ti), 1 ^n, 

Zi,n{Ti-l) = Zj_l,„(Ti_i) + in fl, 2 ^ Z ^ 71, 

^l,n(b) — Zn^n(Y) in H. 


Zi^n — 0 


( 2 . 12 ) 


Since {Y*, U*) is the optimal solution of (lOCP)n, we have Jn{Yn,\, Un,\) — Jn{Y*,U*) > 0. Dividing 
by A and passing the limit A —>■ O'*", using (1.6), (2.11) and (2.12), we infer that 

n pn 1 ^ 

E/ {yln-yd,Z,,n)dt+—Y,{uU,n.U^-l,u)>Q, V C/„ G (L^ (D))”"!. (2.13) 

i=i dTi-1 

Let p* be defined by (1.9) (same reasoning as in Section 2.1). Multiplying the first equation of (2.12) by 
Pn and integrating over D x (Ti_i,Ti), we get 


{Zi^n{Ti),pl{Ti)) - {Zi^n{Ti-l),P*nin_i)) = / (?/*„- J/rf, dt, 1 < Z ^ 71, 

-'Ti-l 
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and summing over i = 1, 2,..., n, we obtain 


n n 

X! / iyln-yd,Zi,n)dt = -'^{Xupl{Ti-i),Ui-i^n), (2.14) 

i=l i=2 

which, combined with (2.13), yields (1.10). 

Let us now prove the converse, that is, let us prove that, if {Y*,U*) and p* satisfy (1.8)-(1.9)-(1.10)- 
(1.11), then {Y*,U*) is the optimal solution of (IOCP)„. 

Given any t7„ = (mi_„, it 2 ,n, ..., Wn-i,n) £ (L^(f2))"“\ we denote by Y„ = (yi,„, y 2 .n,..., 2/n.n) the 
corresponding solution of (1.7). By using arguments similar to those used to establish (2.14), we obtain 
that 

n n n 

/ ^ylr>^ ~ yd^ y^n - vln) dt + ^{XujPniTi-l), Mi-1,n “ = 0. 

i=l i=2 

This, together with (1.10), implies that 


Hence 


E 

2=1 ‘ 


pTi 2 

I {yin - Vd, Vi,n - vln) + 7 — Mi- 1 ,n “ M*_i „) = 0. 

^ 2=2 


JniYn,Un)-Jn{Y:,U:) 

= o E / ^2/i,n + yin - ‘^yd,yi,n - vln) - M*_i „) 


i=l 


2=2 


^ 'y / (2/i,n yd,yi,n 2 /i,n) dt + , Mi-l,7i Mj _2 „) — 0. 

i=l dTi-i n .^2 


We conclude that {Y*,U*) is the optimal solution of (lOCP)n. 


2.3 Proof of Theorem 1 
2.3.1 A first estimate 


The following lemma compares two states generated by controls activated in different ways. 
Lemma 2. Let 0 ^ Ti < Ti + S < T 2 < +00 and let u € L^(fl). Let z and w he the solutions of 


and 


dtZ - Az = ■^X(Ti,Ti+ 5 )Xl,;M 
z = 0 


z{Ti) = 0 


in Hx(Ti,T 2 ), 
on dfl X (Ti,T 2 ), 
in LI 


dtw — Aw = 0 
w = 0 

w{Ti) = Xu,u 


in nx{Ti,T 2 ), 
on aHx(Ti,r 2 ), 
in LI. 


For every p G [2, 00 ), there exists C{T 2 ,p) > 0 such that 


\\z - w\\LP{Ti,T 2 ;L^{n)) ^ Ct{T 2 ,p)Si’ ||XtjM||. 


(2.15) 


(2.16) 
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Proof. By the definitions of 2 and in, we have 


{t)-w{t)= j ^)ix(Ti,Ti+ 5 )('r)Xa;MdT-e^(* Vte[Ti,T 2 ]. (2.17) 


Let q = and let / G L^(Ti, T 2 ; L^(il)). Let (p be the solution of 

dt(p + Aip = f in nx{Ti,T 2 ), 
if = 0 on dn X {Ti,T 2 ), 

^p(T 2 ) =0 in f2. 

By [19, Theorem 1], there exists C{T 2 ,p) > 0 such that 

\\dtp\\L<i(Ti,T2;L^m) ^ C!{T2,p)\\f\\Li{TuT2;L2{n))- 

It follows from (2.17) and (2.18) that 


(2.18) 


(2.19) 


pT2 / 1 /*Ti+(5 

/ {z{t) - w{t)J{t))At = lxu,u,p{Ti) - - j pij) 

JTi \ ^ JTi 


^Ti+(5 pT 


dr ) = - ( XojU, ■ 


dtp dt dr ) , 


/Ti JTi 


which, together with (2.19), yields 
rT2 


iTi 


{z{t) - w{t)J(t)) dt 


^ IIXi.c;U|| [ \\dtp\\dt C{T2,p)6^^^\\XuJu\\\\f\\L1{T^,T2■,L^n))■ 
JT^ 


This leads to the desired result and completes the proof. 

2.3.2 Proof of the control error estimate 

In this section, our objective is to establish (1.19). 

Recalling that u* is defined by (1.18), we denote by y{u’^) and by p(u*) the solutions of 

dtviK) - ^y«) = XujU* in 11 X (0, T), 
y('^n) — d on dfl X (o,r), 

y«ko) = y(u*}(T) in n 

and 

dtPiK) + ^piK) = yiK) - yd in x (o, t), 

P(“n) = 0 on 9Hx(0,T), 

P(<)(0) =P«)(r) in H. 

The existence of these solutions follows from Section 2.1. The proof goes in three steps. 


□ 


( 2 . 20 ) 


( 2 . 21 ) 


Step 1. We claim that 


dt = /i+/2 


( 2 . 22 ) 


with 


^ ^ pTi pTi 

dl = X! / ~ XujP{K),U* - O dt, -^2 = X! / ixujpiul) - XujPn{n-i),u* - <) dt, 
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where p* and p* are given by (1.4) and (1.9) respectively. 

The claim follows from (1.5), (1.18) and (1.10), and from the fact that 

pT pTi ^ pTi 

/ Ik* - ^ / (M*-<,U*-Odt = ^ / {XojP* - Xu:Pn{n-i),U* -Ul)dt. 

J T\ i—^, ^ —1 i—2 ^ ^* — 1 


II ^ C{T)hn\\yd\\'L2(^o^T;L^{Q)) ■ 


Step 2. We claim that 

We first infer from (1.18) that 

n „Ti 

dl = X! / iXu:{P* - PiK)),U* - U'^)dt 


i=2 
rT 


= [ (P* -piUn),Xu:iu* -U^))dt- f {p* -p{u*„),Xu:U*)dt. 

Jo Jo 

Then, on one hand, by (1.3), (1.4), (2.20) and (2.21), we get that 

dt{y* - y{K)) - A{y* - y{u*J) = Xc^iu* - u*J in flx(0,r), 
y*-y{Un)=0 on dXlx{0,T), 

iy*-yiK)m = iy*-y{u*J)iT) in fl 


(2.23) 


(2.24) 


and 


dt{p* - p{K)) + A{p* - p{u*J) = y* - y{ul) in flx(0,T), 
p*— p(m*)=0 on X (0,T), 

(p* - p«))(o) = (p* - piK))(T) in n. 


(2.25) 


(2.26) 


Multiplying the first equation of (2.25) by p* —p(u*) and integrating over x (0, T), by (2.25) and (2.26), 
we obtain that ^ ^ 

[ {p*-p{u*J,xU^*-K))dt = - [ ||j/*-2/«)f dt<0. (2.27) 

Jo Jo 

On the other hand, since {Y*, U*) is the optimal pair for the problem (lOCP)n, we have Jn{Y*, U*) ^ 
J„(0,0). Then by (1.6), (1.12) and (1.18), it follows that 

r Mr dt+ T||<f dt^C Tlly.f dt. (2.28) 

Jo Jo Jo 

From (2.28), (2.20), (2.21) and Lemma 1, we infer that 

l|j/(^n)llc([ 0 ,T];ffi(n)) + Ib(Wn) II C([ 0 .T] i (n))nff i ( 0 T;L 2 (n)) ^ C'(T’) lll/dll L 2 ( 0 ,T;L 2 (n)) . (2.29) 

Since {y*,u*) is the optimal pair for the problem (OCP), we have J(y*,u*) ^ J(0, 0). Then by (1.1), 
we get 

T|b*fdt+ r WuMdt^C r WyMdt, 

Jo Jo Jo 

which, combined with (1.4) and Lemma 1, implies that 


lb*llc([ 0 .T];L 2 (n)) ^ C'(T)||yd||L2(o,T;L2(n))- 


(2.30) 

(2.31) 
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By (2.24), (2.27), (1.5), (2.29) and (2.31), we get that 

h ^ - f {p* -p{ul),Xu>P*)<it 

Jo 

^ nib* -p(u*)||c([0.T];L2(n))lb*llc([0.T];L2(a)) < C'(7’)/ln |bd|li2(o_T;L2(n)) > 

and (2.23) follows. 


Step 3. We claim that 


Ii ^ 


■ f lb* — dt + C'(T)/l„|b(i|||2(o^T;L2(n))- 

J n 


We first note that 


-^2 = X! / (x<.^7’«)-X<^Pn(n-i),n*-Odt 

i = 2 




-<f dt + i^ / lb«)-K(n-i)f dt. 


i=2 


(2.32) 


(2.33) 


Then, we proceed with three sub-steps. 

• Sub-step 3.1. Let us prove that 

^ PTi pT 

E/ lb(<)-p:(n-i)f dt^C'(r)/r2||y,||2^(^_^,^,(^^^+C(T) / |b«)- 2 /:fdt. (2.34) 

i=2 -^0 

By (2.29), we have 


E/ lb«)-P«)(n-i)f dt = E / / bp«)ds dt 

2—2 — 1 2=2 — 1 ‘^i — 1 

^ ^nllbp(w„)|b2(g ^ C'(T)ft,„||j/d|b2(o,T;L2(n))- (2.35) 

Moreover, from (1.9), (2.21) and Lemma 1, we obtain that 


E / lb«)(n-i)-p);(n-i)f dt T|b«)-p);||c([0.T];L2(n)) ^ C'(T) / |b«)-y;fdL 

i=2 •'’"i-i •'0 

Combined with (2.35), this gives (2.34). 

• Sub-step 3.2. We claim that 

l|e‘^(2/(w)l)(0) - 2/^(0))||LP(o,T;L2(n)) < C(T,p)hl/P\\yd\\L^(^o^T-,L^{Q)), (2.36) 

for every p G [2, oo). 
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Indeed, by (1.8), (1-11), (2.20) and (1.18), we have 


2/»(0) = e'^^J/n(0) + 


J=2 


y«)(0) = e^'^y«)(0) +J2ir di- 


i=2 


Then 

J/«)(0) - 2/:(0) = e^^(y«)(0) - 2/:(0)) 

' 1 




,=2 \ " 




(2.37) 


It follows that 
tA 


I|e (2 /(m*)( 0) - ?/*(0))||z,P(o,T;L2(n)) 


^ (l_e-^i^)-i^ 

i=2 


hn Jr,., 


LP(0,T;L2(0)) 


^ c(r)^ 

i=2 


'-n Jt,_i 


p(‘-'^)^v ?;* dr - pd-'>-i-i)A * 


LP(rj,2T;L2(n)) 


where —Ai < 0 is the first eigenvalue of the Dirichlet Laplacian. Moreover, by Lemma 2 with 
Ti = Tj-i, 5 = hn, T 2 = 2T and u = u*_^ we get that 


/ e(‘ ^)^X(r,_i.r,-)(T)Xc^Uj-l,ndT- ed 

JTj-i 


LP(rj_i,2T;L2(n)) 

< C{T,p)hl/P\\xu,u*-iJ, (2.38) 


for every j G {2,..., n}. Since 

f e(*“^)^X(r,_i,r,)('r)x<..wJ_i,„dT = / ed-^)^Xc^M*_i^„dT, VtG[r,-,2T], 

Jrj_i 

by (2.38) and (1.18), we obtain that 

n 

||e‘^(y«)(0) - y:(0))|Up(o,T;L^(n)) < C{T,p)Y,hl!^uU,n\\ 

1=2 

^ C{T,p)hl!^ I ('Zh\uU,n 


(2.39) 


^1=2 


Vl=2 


^ C{T,p)h]/^ KllL2(0.T;L2(n)), 
which, combined with (2.28), implies (2.36). 
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(2.40) 


• Sub-step 3.3. Let us prove that 

l|j/(K) - ynl|L!>(0,T;L2(n)) ^ (^(T,||T/d||L2 (o_T;L2(O)), 

for every p € [2, + 00 ). 

Let (zj^n)i^j^n and (wj^n)i^j^n be solutions of (2.15) and (2.16) respectively, with Ti = 
hn, T 2 = T and u = We set 




^j,n(t) — 


t e (o,T,_i], 


(L? — 

By (1.8), (1.11), (1-18) and (2.20), we have 

y«)W - Vnit) = e‘'^(y«)(o) - y)^(o)) 


and Wj^n(t) = 


t e (o,T,_i], 

j,n{t), t € (rj_i,T). 


+ X(r,_i,t)(s)e(* '’)^x„<(s)ds-e(* 

j=i Y^j-i 

i 

= e‘'^(2/«)(0) - yn(0)) + Y^{Zj,n{t) - 

i=i 


(2.41) 


for every t G (r^-i, r^] and every i G {1, ..., n}. Then 

n 

y«)W - Vnit) = e‘'^(2/«)(0) - 2/0(0)) + '^XY_,,T){t) {zpn{t) - Wj,u{t )), (2.42) 

1=1 

for every t G (0, T). Indeed, given any t G (0, T), let io G {1,..., n} be such that t G (r^o-i, Ti^]. It 
follows from (2.41) that 


2/«)W - yOW = e*^(y«)(o) - y 0 (o)) + 

1=1 

io 

= e*^(y«)(o) - y0(o)) + 


1=1 

^0 

= e*^(y«)(o) - y0(o)) + iAn{t) - A At)) 

1=1 

which yields (2.42). By (2.42), (2.36) and Lemma 2, we obtain that 

lly(llO) - 2/0 IUp(O,T;L 2(n)) 

n 

^ ||e‘^(j/(uO)(0) - 2/0(O))||LP(O,T;L2(a)) + '^Aj,n - Wl.n||LP(rj_i.T;L2(n)) 

1=1 

^ C{T,p) ^n^\\yd\\L^{0,T-,L‘^{^)) + j2f^A\\nAiA\^ ■ 

Using (2.28) and the same arguments as in (2.39), we obtain (2.40). 

Step 3 follows immediately from (2.33), (2.34) and (2.40). 

Finally, the theorem follows from (2.22), (2.23), (2.32), (1.18), (1.5) and (2.31). 


19 


2.3.3 Proof of the state and cost functional error estimates 

In this section, our objective is to establish (1.20) and (1-21). 

We start with the case 2 ^ p < +c». By the triangular inequality, we have 

\\y* - yn\\LP{0,T-,L^{n)) ^ \\y* - yiUn)\\LP{0,T-,L^in)) + ||y«) - yn\\LP{0,T-,L‘^(n))- (2.43) 

We infer from (1.3), (2.20) and from Lemma 1 that 

\\y* - ?/«)IIl.(0.T;L=(O)) < T-p \\y* - y«)||c([0.T];L=(a)) < C{T,p)\\u* - <|U2(0,T;L^(O)). (2.44) 

Since p ^ 2, it follows from (2.44) and from the first part of Theorem 1 that 

lly* “ 2/(M^)l|LP(o,T;L2(n)) ^ C'('r,p)\/^||yd||L2(o.T;L2(n)) ^ C!{T,p)hl/^\\yd\\L^(o,T-,L^in)), 

which, combined with (2.43) and (2.40), gives (1.21). 

Finally, (1.20) follows from (1.1), (1.6), (1.12), (1.18), (1.19), (1.21), (2.28) and (2.30). 


2.4 Proof of Theorem 2 


2.4.1 A general result in measure theory 

Lemma 3. Let io be a measurable subset of fl having a boundary dtu. For e > 0, we define 

uje = {x & I d{x, duj) < e}, (2.45) 

where d{x,duj) = inf{|a; — y\ \ y € duj}. There exists p > 0 such that 

\dujrj\dr] ^ 2{1 + e/p)^~^\duj\e, (2.46) 

for every e € (0, p). 

In (2.46), without ambiguity, [weI designates the Lebesgue measure of uje, and \dujn\ = 'H^~^{duJr,) 
designates the {N — l)-Hausdorff measure of dujn- 

We give a proof of this result for completeness, borrowing arguments from [11]. 

Remark 6. In the proof below, the assumption duj G is required. For the general case, whether 
(2.46) holds or not seems to be open. 

Proof. For every y G dtv, let v{y) and r(i/) respectively denote the unit inner normal to dtv at y and 
the tangent hyperplane to duj at y. The curvatures of duj at a fixed point j/o G du are determined as 
follows. By a rotation of coordinates, we assume that the xn coordinate axis lies in the direction v^yo). 
In some neighborhood Af{yo) of y^, we have MIjjq) H doj = {xn = ^{x')}, where x' = (a;i,..., xn-i), P G 
C'^(r(j/o) nA/’(?/o)) and D(p{yQ) — 0. The eigenvalues ki, ..., kn-i of the Hessian matrix Z1 ^(/j(j/q) are the 
principal curvatures of dtv at yo and the corresponding eigenvectors are the principal directions of duj at 
2 / 0 . By an additional rotation of coordinates, we assume that the xi,... ,xn-i axes lie along principal 
directions corresponding to /ci,..., kat-i at yo. Such a coordinate system is said to be a principal 
coordinate system at yo. The Hessian matrix with respect to the principal coordinate system 

at yo described above is given by D'^p(yQ) = diag(/ti,..., kn-i)- The unit inner normal vector v(y) at 
the point y = (y', p{y')) G Al'(yo) H div is given by 



^(y) = - 


DzPiy') 


y/lHWF’ 


1 ^ ^ A- 1, 


J^w(y) = 


^l + \Dp{y'W' 
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Hence, with respect to the principal coordinate system at j/Qj we have 

DjVi{yo) =-KiSij, j = 1, ■ • ■ ,-/V - 1. (2.47) 

Since dui is dui satisfies a uniform interior and exterior sphere condition, i.e., at each point yo G dui, 
there exist two balls Bi and B 2 depending on yo such that Bi fl — uj) = {j/o} and i ?2 fl ZU = {j/o}) 
and the radii of the balls Bi and B 2 are bounded below by a positive constant denoted by y.. It is easy 
to show that bounds the principal curvatures of duj. 

The rest of the proof goes in two steps. 

Step 1. Let us prove that {Q < e < y) has a C^-smooth manifold structure. 

Given any point x such that d{x,duj) < y, there exists a unique point y = y{x) G dui satisfying 
\x — y\ = d{x, dui). We have x = y + v(y)d{x, duj) if a; G w and x = y — u(y)d{x, dui) ii x ^ u. Now we 
give a construction of a C^-smooth manifold structure on uig. For this purpose, we fix a yo G du and we 
define the map $0 from U = (r(?/o) nA/’(yo)) x (—to by 

‘^o{.y\d) = y + iy{y)d, V {y',d) € {T{yo) nJ\f{yo)) x {-y, y), (2.48) 

where y = (y', ^{y'))- By (2.47), the Jacobian matrix of 4>o at (yg, d) is I?$o(yo, d) = diag(l —Kid,..., 1 — 
KN-id, 1), and hence det iJ$o(yo, d) = (1 — Kid) • • • (1 — Kjv-id) 7 ^ 0, for every d G {—y,y). It follows 
from the inverse function theorem that $0 is a local C^-diffeomorphism in a neighborhood of any point 
of the line {yp} x {—y,y). Then by compactness of [—e,e], we can choose Uq = Bq x [—£,£], with Bq 
an open ball in r(yo) fl A/’(yo), so that $0 is a C^-diffeomorphism from Uo to $o(ZJo)- This shows that 
($o(ZJo), ^) is a coordinate chart centered at yo in the topological space Wg. 

We carry on the above process for each y G dui and we define an atlas {(14, 'h^^)} for Wg, where 14 is 
an open neighborhood of y^ G dui, ^(14) =Ua = B^x [—e, e] and Ba is an open ball in r(y(j) r\M{ya). 
By the definition of (similar to (2.48)), one can check that any two charts in {(I4,<I>“^)} are C^- 
smoothly compatible one with each other. Hence {(14, is a atlas for uig. This atlas induces a 

structure on Wg. 

Step 2 . Let us establish (2.46). 

By [II, Lemma 14.16, page 355], we have d(-, 9a;) G C^(a;e) and |Vd(-, 9a;)| = 1 in uje, which, combined 
with the Coarea Formula (see, e.g., [9]) applied to / = d{-,dui), gives 

|Vd(a:,9a;)| da: = [ n^~^{{d{-,duj) = rj}) drj = f |9a;,,| dy = f (|9a;+| + |9a;“|) dy, (2.49) 
Jo Jo Jo 

where dui'l and duij are the inner and outer parts (with respect to a;) of duis for each S G (0,6:), 
respectively. Now, given any y G (0,e), to compute |9a;+| and |9a;“|, we define ipr) ■ duj ^ 9a;+ by 
iirjiVa) = ° Trj o $“^(ya), for every ya G duj, where is the mapping given by t^(z,0) = {z,ii) for 

every z G From Step 1, we take two arbitrary coordinate charts {{Vp,^~^^)} and {( 14 ,$“^)}, 

where 14 and 14 are open neighborhood of yp and y.^ ( y/3i 2/7 £ respectively. Then by the definitions 
of and (similar to (2.48)), one can check that 

^/3 o o $“4y) = ^7 o 4 o Vy G Vs n foy n9a;. (2.50) 

We recall from Step 1 that each is diffeomorphic from I 4 to Ua = ^a{Va). Therefore, by (2.50), 
iprj is diffeomorphic from duj onto dui^ and 

det($„ oTrjO $”^)(ya) = (1 - Ki{ya)r]) • • • (1 - KN_i{ya)r]) G ((1 - (1 + e/r"^)^"^), 
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for every ya £ duj. This, together with the definition of 'ipri and (2.50), implies that 

= / det['(/',,](a:) dtr ^ (1 + I, ^ y £ {0,e). (2-51) 

J duj 

Similarly, we have \duj~ \ ^ (1 + r]yL~^)^~^\du:\^ for every rj £ (0,e). Then, (2.46) follows from the latter 
inequality, from (2.51) and (2.49). 

This completes the proof. □ 

2.4.2 Smooth regularizations of characteristic functions 

We define the C°° function xh ■ —>■ K by 


for every x £ where 


Xlix) 



y)xuj{y) dj/, 





ijix) 


if |a;| < 1, 
0 if |x| ^ 1, 


with c > 0 such that r]{x) dx = 1. 

Lemma 4. Let y he as in Lemma 3 and let e £ (0, /i). For every p £ [1, +oo], we have 
llVx^llLP(n) ^ Ce~^^p and \\xl - xAlp(xi) < Cep. 


(2.52) 


(2.53) 


(2.54) 


Here and throughout the proof, C is a generic positive constant independent of p and e. 

Proof. Note that the case p = +oo follows by passing to the limit. Therefore it suffices to prove (2.54) 
for 1 ^ p < + 00 . We set uj^ = {x £ uj : d{x, doj) > e} and = {x ^ to ■. d{x, doj) > e}. Then ujI and 
are open subsets of such that 

w] Ua;2 Uwe = R^, (2.55) 

where Wg is defined by (2.45). 

On the one hand, by (2.52) and (2.53), we get that 

xUx)=f -^v(- — -')xuj{y)dy, (2-56) 

Js.n £ \ £ J 

which, combined with (2.53), yields 


'^xlix) = [ ^-Xa;(y)Vr7 


x-y 


dy = - / Xu;{x- £y)Vp{y) dy 

£ JR'V 


-f 

£ J{yem^-.IvKi} 


Xc,;(x-ey)V? 7 (y)dy. (2.57) 


On the other hand, by (2.56), we have 

Xlix)= r]{y)xu,ix - ey) dy = 

Jrw J( 


{yeR":|y|<l} 


y{y)Xi.^{x - ey) dy. 
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This implies that 


Xl{x) = 


1 if X G Wg, 
0 if X G Lo'i. 


(2.58) 


It follows from (2.55), (2.57) and (2.58) that ||Vx^||^p(q) ^ |Vx^(x)|*’dx < {Ce ^)p|w£|, which, 

combined with Lemma 3, yields ||Vx® ||iP(o) ^ Ce~^^p . 

Besides, by (2.58) and (2.55), we have 


\xUx) - xUx)\ < 


0 if X G Wg U Wg, 
1 if X € ujs- 


This, together with Lemma 3, implies that \\xli~Xc^\\Lp{n) ^ \^s\^ ^ Cep . This completes the proof. □ 

2.4.3 A useful estimate 

The following estimate for a linear heat equation is not standard. 

Lemma 5. Let lu C LI be a subset having a boundary. Let p G (l,+oo), let Ti and T 2 be two 
nonnegative real numbers such that Ti < Tj, and let zq G Hq{LI). Let z be the solution of 


dtz — Az = 0 in LI x {Ti,T 2 ), 
z = 0 on dLlx{Ti,T 2 ), 

z{Ti) = XujZo in fl. 


\\z{s)-z(Ti)\\ < (s-Ti)2||zo||ffi(o), 


If uj = LI, then 

for every s G [Ti,T2]. 

If uj LI, then 

( C{T2){s-T,)^\\zo\\H^,in) */ 

||z(s)-z(Ti)|K C{T2,p){s-T^)^p\\zo\\Hiin) */ N = 2, 

[ C{T 2 )is-T,)i\\zo\\Hiin) N = l, 

for every s G [Ti,r2], for some constants 0X2) > 0 and C(T2,p) not depending on zg. 


(2.59) 


(2.60) 


(2.61) 


Proof. Since the proof of (2.60) is similar to obtain but simpler than the one of (2.61), we assume that 
we are in the (more difficult) case where w LI. Let p be as in Lemma 3 and let s G (Ti, T 2 ]. We set 


Co = 


min{p,, 1} 


and e = cq v^s — Ti. 


2 max{i/72, 1 } 

Note that e < min{^, 1}. By (2.59), we have 

z{s) = e(^-^^)^(x. - X^)^o + = zi(s) + ^2(5), 

for every s G [Ti, T 2 ]. We have 


(2.62) 


(2.63) 


\\Z2{S) - Z2in)r = 


dtZ2 dt 


Ti 


dx^{s-Ti) \\dtZ2\\‘^dt, 

J Ti 
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for every s G \Ti , T 2 ]. By definition, Z 2 is the nnique solution of the Dirichlet heat equation with initial 
condition z^iTi) = By integration by parts, we have 

[ \\dtZ 2 \\'^dt= f f {dtZ 2 {t,x))'^dxdt= f [ dtZ 2 {t,x) ■ Az 2 {t,x) dxdt 

JTi JTi Ja JTi JQ 

= i f \Vz2{T„x)\^dx-^ j |Vz2(T2,x)|2da:^ ||V22(Ti)f, 

\J O t/ o 

and therefore we get that 

||^2(s)-Z2(ri)f <(s-ri)||VZ2(Ti)f, 
for every s G [TijTlj]. It follows from (2.63) that 

\\z{s) - z(Ti)|| ^ ||zi(s) - z,{T,)\\ + ||Z 2 (S) - Z 2 {T,)\\ < 2||(x^ - xl)zo\\ + V(x^zo)||. (2.64) 

2N 

N ^ 3, then, using the Holder inequality, the Sobolev embedding ^ Lemma 4 and 

(2.62), we obtain that 

IKXo; - xt)zo\\ < - xlWh^in) < C'||^o||i?i(n)e^, 

and 

IV(Xw^o)|| ^ Ikollffi(n) + ||^:o||^,2^^j.^^||Vx5llLW(f2) ^ C||zo||//i(o)(l + £« ^). 

These estimates, together with (2.64), imply that 

||z(s) - z(Ti)|| < C\\zo\\Hiin) + Vs-Tie^-i). (2.65) 

From (2.65) and (2.62) it follows that (2.61) holds. 

If = 2, then, similarly, 

IKXo; - Xl)zo\\ ^ Ikoll^j^^^^llXo; - XlWL^Pin) < C'(p) 11^0 ||ff i (n)^^, 

IV(Xw^o)|| ^ lko||i?i(n) + II^oI1^^^^^||VxS,||l2p(q) < C'(p)ll^o||ffi(n)(l + e'^p ^), 
and using (2.64) we infer that 

\\z{s) - z{Ti)\\ < C'(p)||zo||ffi(n)(e^ + \/s - (2.66) 

It follows from (2.66) and (2.62) that (2.61) holds. 

If = 1, then 

IKXc. -X^)^o|| < lko|lc(n)llx<^ -xLII < C'lko|ki(n)e^ 

IV(Xw^o)|| ^ INo||_H-i(n) + ll^o|lc(a)IVx!;ll ^ C'lko||//i(n)(l + e 

which, combined with (2.64), imply that 

|| 0 (s) - z(Ti)\\ < C||zo||//i(n) _|_ y'g _ Tie =). (2.67) 

By (2.67) and (2.62), we obtain (2.61). The proof is complete. □ 
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2.4.4 Proof of the state error estimates 

We prove (1-22) only when N ^ 3, the other cases being similar. Let u* and U* be the optimal controls 
solutions of (OCP) and (lOCP)n respectively, where U* = ^ 2 .Let 

M* be given by (1.18). We have 

\\y* - y))l|L“( 0 ,T;L 2 (a)) ^ \\y* - y(u*)||Loo( 0 _T;L 2 ( 0 )) + \\y{ul) - 2/*||L~(0,T;L2(n))- (2.68) 

By (1.3), (2.20) and Lemma 1, we infer that 

\\y* - y{Un)\\c{[0,T]-L^(n)) ^ C{T)\\u* - Un\\L^{o^T-L'^{Q.))- (2.69) 

Besides, we claim that 

WyiKi) - 2/nllL“(0,T;L2(n)) < \\yd\\L^{ia^T-L'^{a))- (2.70) 

Then (1.22) follows from (2.68), (2.69), (2.70) and Theorem 1. 


Let us prove (2.70). On the one hand, by (2.37) and (1.10), we have 

1 , ^ 


ll2/:(0)-y«)(0)K 




1=2 


h 


n 


(T-s)A * Ag_JT-Tj-i)A 


^C[T) 


n pTj 

^ (/ - Xu.p:{t,-i) ds 


(2.71) 


< C'(T)^ / (j - x.p:(r,-i) 

1=2 ^ ' 


1 = 2 -' O - 

On the other hand, by (2.41), (1.18) and (1.10), we infer that: 
• for every t G [0,ti], 


ds. 


(2.72) 


2/«)W - vlk^) = e*^(y«)(o) - 2/i^(o)); 

for every t G ('ri,T 2 ], 

y{u*n){i) - y*u{t) = e*'^(2/«)(0)-y);(0))+ / (2.73) 

'J T^ 


' Ti 

rt 


= e‘'^(l/«)(0) - yi)(0)) + / '*)^x„p*(Ti)ds- f ds; 

J Ti J Ti 

for every t G (ri_i,Ti], with 1^3, 


y{<)it)-yUt) = e*^(y«)(0 )-yn(0))+^(/ '*)^x^<(s)ds-e(‘ 

j =2 

+ f e(‘"")^x^<(s) ds - (2.74) 

Ti-l 


2-1 


e‘^(y«)(0) - y:m + ^ - e(^---^)^]x.p:(r,_i) ds 

+ [ e^*~''^'^XcjP*nin-i)ds- [ ed-'^'-d^^^pl(Ti-i)ds. 
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It follows from (1.9), Lemma 1 and (2.28) that 

lbnllc([0,T];//i(n)) ^ C")?") II J/d || (0,T;L2 (Q)) , 

which, combined with (2.72), (2.73), (2.74) and (2.71), implies that 


ll2/(Wn) - 2/nl|L“>(0,T;L2(n)) 

pT. 




r'^3 

C{T)Y, / II [/ - e(^-"-^)^]x.p:(r,-i)|| ds + C(T)h„||2/,|U.(o,T;L^(a)). 

- 7=9 'Tj-l 


j=2 ^3 

This, together with Lemma 5 and (2.75), yields 
\\y{Ki) - 2/nllL“(0,T;L2(n)) 


" /•Tj ^ 

^ C'(T)^ / (s - Tj_i)^||p*(Tj_i)||^l(n) ds + C(T)/l„||?/d||L 2 ( 0 ,T;L 2 (O)) 

1=2 "'■j-l 


3=2 ■' 

^ C'(T)hy^^||7/d||L2(0_T;L2(n)) 

and (2.70) follows. This ends the proof. 


(2.75) 


2.5 Proof of Theorem 3 

2.5.1 Proof of the control error estimate 

Let us establish (1.23). As in Section 2.3.2, the proof goes in three steps. 


Step 1. We claim that 


with 


i=i “''^■-1 


iw*-<„f dt = J1+/2, 


n „ri 
i=l “'^>-1 


and 


-^2 = (xujPn - j pl{s)ds,u* - vl^j dt, 

where p* is defined by (1.4) and p* is defined by (1.15). 

The claim follows from (1.5), (1.16) and from the fact that 


(2.76) 


E 


l«*-<nf 




= E/ ~ J Pn(p)ds,u*-vl^j dt. 
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Step 2. We claim that 


(2.77) 


/i <0. 

Indeed, using (1.3), (1.4), (1.14) and (1.15), we get that 


and 


dt{y* - Vn) - - Vn) = Xu,{u* - f*) in rjx(0,r), 

y*-yn = ^ on 511x(0,T), 

{y* - yu){0) = {y* - y*n){T) in 

dt{p* - Pn) + ^{p* - Pi) = y* - Vn in llx( 0 ,r), 
P*~Pl = ^ on dXt X {Q,T), 

{p* - Pl)i0) = {p* - Pl)iT) in a 


(2.78) 


(2.79) 


Multiplying the first equation of (2.78) by p* — pi and integrating over x (0,T), by (2.78) and (2.79), 
we obtain that ^ ^ 

[ {P* - Pl^Xujiu* - fl))dt = - [ ||y* - dt ^ 0, 

Jo Jo 

which, combined with (1.17), gives (2.77). 

Step 3. We claim that 


I/ 2 I ^ C{T)h4ydhHo.T;mn)) J 


Indeed, on one hand, we first note that 

n 

I'-A « E / 

i=l 


Pn - 


plis)ds 


Jn-i 


\u* - <n|| dt. 


(2.80) 


(2.81) 


It is easy to check that, for every t £ r^]. 


1 r' 


1 r' 

1 

r' /■* 

Pni^)-ir Plis)ds 

JTi-i 

— 

T- (Plit) - pUs)) ds 

Jn-i 

II 

^11 
s 1 

/ / '9rP„(T)drds 

J Ti — I J S 




/ T \ 

[ \\9rPliT)\\dTds= f \\drPliT)\\dT i f || (t) f dr ) . (2.82) 

Jn-i Jn-i Jri-I \JTi-i j 


It follows from (2.81), (2.82) and from the Holder inequality that 


I/ 2 I < WdtplW^dtj J ||u*-w;„||dt 




u-vlJ^dt 


< 




(2.83) 
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On the other hand, since {y^,V*) is optimal (with 
from which it follows that 



)), we have J{yUX) ^ 


WVn -ydfdt / Wvdfdt, 


(2.84) 


0 


and 



(2.85) 


By (1.15), (2.84) and Lemma 1, we get that ||9tp*||^dt ^ C'(T) fj' || 2 /d||^dt . This, combined with 
(2.83), implies (2.80). 

Finally, (1.23) follows from (2.76), (2.77) and (2.80). 

2.5.2 Proof of the state and cost functional error estimates 

We start with establishing (1.24). Using (1.3) and (1-14), we have 



0 X (o,r), 
dn X (0,T), 

n. 


( 2 . 86 ) 


(1.24) follows from (2.86), Lemma 1 and (1.23). 

Finally, by (1.1), (1-23), (1.24), (2.84) and (2.85), we obtain (1.25). 
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